﻿ 4 Approximation of continuous time systems with discrete time systems A Theory (it helps if you read the course too) The continuous-time systems are replaced by discrete-time systems even for the processing of continuous-time signals 1) Impulse Invariance Method: () [] dahn ThnT=⋅ The relation between the Z transform and the Laplace transform: ∞ 2k⎛⎞π () sTHz Hsj=− ∑ da⎜⎟ ze= ⎝⎠ ∞ kT=− The frequency response of the digital system is the same with the frequency response of the analog system of limited band for frequency less than half of sampling frequency ππ HH=Ω ≤ ≤; and ωωω () () adT MΩ=ω TT 1 The relation between the s and z planes Tσ ⎧ = jΩsT j T j TreΩσ ω ; s σ = =+ = ⇒ = ⋅ ⇔ jω zre xjy ; ze re e e=+ ⎨ 2Tkω =Ω+ π ⎩ 2) Step invariance method: [sn snT= () ] da 3) The method of the differential equation approximation Differential eq: kk NM dyt dxt ()() ab= kkkk∑∑ k0 0kdt dt== Use approximation: 1yn yndy t y nT y nT T−− ()()()[][] ddaaa−− ≅= tnT= dt T T Finite difference eq: Nk Mk 11pp pp () −− 11aCynpbCxnp−− () [][] ∑∑ ∑∑ kk kkkk= k0 0 0 0pkpTT== == The relation between the Z transform and the Laplace transform: = () () 1daHz Hs− 1z− s= T The relation between the s and z planes 1− 1z− s= T 2 4) Bilinear Transformation Method 1− 1Tz+ () Hz Hs== () 11da− 21z−− 21sz=− 1T− 1z+ B Problems Problem 1 We have a linear time-invariant system with the transfer function: 2s+ Hs= () 2a + 13ss+ ()( ) a) Determine its impulse response, ( ) aht ; b) Find the impulse response of the digital system equivalent to the analog one, using the impulse response invariance method [ hn; ] d c) Find ()Hz, the Z transform of [hn ; d]d 3 Solution a) To determine the impulse response from the transfer function, we decompose it into simple fractions: ABC + () 2aHs=+ 131sss++ + () 221 1s+ = =⇒= 1Bs Hs B=+ () () − 11ass=− + 32 2s = 21s+ = =− ()() 23aCs Hs=+ 33ss=− 41s =− + () −− 22321ddsss++ ⎡⎤ ⎡⎤ = = = () () 21aAsHs=+ ⎢⎥ ⎣⎦ − =− 11 1ss s=− + 343ds ds ss = ⎣⎦ + () We identify each term using the tables of Laplace transforms: 11 1 1 11 +⋅ −⋅ () 2aHs=⋅ 412 431sss++ + () 77 311tt−− et e t−σσ ()() 44 11d⎛⎞ ⎜⎟ 2−= 11ds ss+ ⎝⎠ + () () st stdX s∞∞ −− = − =− =− xtXs xttedt txte txts↔⇒L () ( )()( ) () ()() {} ∫∫ −∞ ds−∞ 11 tt−− − ⇒−↔−σσL ()()() {} 22te t s te t−= 11ss++ () () 311 1ttt−−− −σσσ ()()()() aht e t te t e t=+ 42 4 2 3nT nT nTTT T−− − b) []=⋅ = + −σσσ () [][][] dah n T h nT e n ne n e n 42 4 c) We identify each term using the tables of Z transforms: /4nTT T nT T −− =↔σσ [][] () 1Ten e n−− 44 1ez− 33/4nnT TTT T−− − ↔−σσ en e n−= [][] () − 31T− 44 1ez− 22 21T−− nTT Tez nT T−− =↔σσ [][] () 12Tne n n e n−− 22 2(1)ez− 4 21T−− /4 /4TTe z T ()=+⋅ − 11231dTTTHz−− −− − − − 12(1)1ez ez e z−− Problem 2 We have () the transfer function of an analog system, linear and time- aHs invariant, 1 Hs= () 2a 12ss++ ()( ) a) Determine its impulse response, ( ) aht ; b) Compute the impulse response for the digital system in the case of the impulse invariance method c) The step response is () hn ] ast For the digital system, its impulse response is [ d and its step response is [ sn Assume we used the impulse response invariance ] d method to approximate the analog system with a digital one: , () [] adhn ThnT=⋅ Verify that the relation: n () [] ∑ adsn hkT= k=− ∞ is true d) Compute the step response of the digital system, using the step response invariance method: () [] adsn snT= Solution a) To determine the impulse response from the transfer function, we decompose it into simple fractions and we identify each term using the tables of Laplace transforms: 1ABC ++ () 22aHs== 1212 2ssss s++ ++ + () ()( ) 1 = = ()() 211aAs Hs=+ =− 11ss=− 2s+ () 21 = =− 21Bs Hs=+ ()() 22ass=− + 1s =− 211dd⎛⎞ ⎡⎤ = =− =− ()() ⎜⎟ 221aCsHs=+ ⎣⎦ =− 222sss=− + 11ds ds ss =− ⎝⎠ + () 22ttt−−− => () −σσσ ()()() aht e t te t e t=− b) Impulse invariance method The impulse response of the equivalent digital system is: 5 22 2nT nT nT−−− = − −σσσ () [][][][] dahn Thnt Te n nTe n Te n=⋅ c) Impulse invariance method The step response in discrete-time is the response to the unit step: nn∞ = −= =σσ () [] [] [] [][ ] [] ∑∑∑ addd dsn nhn hk nk hk T hkT=∗ kkk=− ∞ =−∞ =−∞ ∞ Which means that: [] The relation is true in the general case () ∑ adsn T hkT= ∞ k=− d) Step invariance method Let ()Ss be the Laplace transform of the step response for the analog system: a () 1aHsAB C D =++ + () 2aSs== ++ + 12 1 22ssss ss ss+ ()( ) + () We compute the step response of the analog system using the decomposition in simple fractions and tables of Laplace transforms: 11 AsSs== = () 00ass== 12 4ss++ ()( ) 1 = =− ()() 211aBs Ss=+ =− 11ss=− 2ss+ () 211 = = 2Cs Ss=+ ()() 22ass=− − + 12ss = () 21213dds+− ⎡⎤ =+ = =− =−= 23/4DsSs⎡⎤ ()() 22a⎢⎥ ⎣⎦ 222sss=− =− + 4ds ds s sss =− + ⎣⎦ () 22113ttt−− − =>() + +σσ σ σ () ()()() ast tetetet=− 424 22113nT nT nT−− − − + +σσσσ () [][] [][][] dasnsnT nenenen== 424 Problem 3 We have an analog, linear time-invariant system with the impulse response hn equivalent to the analog system () ] aht and a digital linear time-invariant system [ d with the impulse response invariance method t− find the transfer function of the analog system; a) If () σ () aht e t= b) Find the transfer function of the digital system; c) Sketch the amplitude-frequency characteristic of the analog system and the discrete-time system Solution 6 1 t− a) ()Hs e t==σL () {} a 1s+ nTT− b) [] hn ThnT Te n H z=⋅ = ⇒ =σ ()() [] 1ddTa−− 1ez− T HΩ= c) () djT−Ω − 1ee−⋅ 11 d) ()HH=⇒ =ωω () aa 21j+ ω ω 1+ TT HΩ= = () 2dTT−− −⋅Ω+ ⋅Ω 221 cos sinTTeje−− 1 cos sinee− ⋅Ω+ ⋅ Ω () T = 2TT−− 12cosee+− Ω For T=0 1 the amplitude-characteristic are almost the same for frequency less than half of sampling frequency 2π/T=10π=31 14 H(ω) 1a H(Ω)|Ω=ωT d 0 8 T= 0 1s 0 6 0 4 0 2 0 -30-20-100102030 Plot it using Matlab What do you notice when the frequency is higher than half of the sampling frequency? 7 >> w=-200*pi:0 01*pi:200*pi; >> Ha=1 /sqrt(1+w ^2); >> Hd=0 1 /sqrt(1+exp(-0 2)-2*exp(-0 1) *cos(0 1*w)); >> plot(w,Ha,w,Hd); grid on Problem 4 Consider the analog, linear time-invariant system with the differential equation: dy 2ytxt+= () () dt a) Find the impulse response and the transfer function of a discrete-time system that approximates the analog system using the impulse invariance method; b) Find again the impulse response and the transfer function if we use the differential equation approximation; c) The same for the bilinear transformation method d) Give implementation forms for the systems from a), b) and c) ; sinxtt= and a sampling period of 0,1T= find the amplitude of the output e) For () signal for the analog system and for the digital system, for the three cases, respectively Solution a) impulse invariance method () 1Ysdy + = ⇒ == 2 yt xt sY s Y s X s H s+= 2 () () () () () () a⇒ 2dt X s s+ () 22tnTT−− σσ = ⇒= ht e t hn Te n H z=⇒ () ()() [] [] 21addnT−− 1ez− b) differential equation approximation T 11 + 112THz Hs− ==== () () − 111zads− − = − +− 11112TzTz− ⋅ 12z− + 12TT+ n 1T⎛⎞ The impulse response of the digital system is: []σ hn n= [] d⎜⎟ 12 12TT++ ⎝⎠ c) bilinear transformation method 1− = ()() 21zaHz Hs− ds= 1T− 1z+ 1− 111z+ === () 11dHz− 11−− 21212zzs−− − + 21 1szT z= −+ + () () + − + 1112Tz− 1Tz+ 8 11−− 111zz++ Hz==⋅ () − 11d− 21 21 2 1 1TTz TTz+− − +−− () ()() 11 1− =+⋅ 1− 21 1 1 1TTTTz−−+−− ()() 11 1 =− + ⋅ − −−+ 1121 1 1TTTT− ()()() 1z− 1T+ The impulse response of the digital system using the bilinear transformation: n 111T− ⎛⎞ δσ hn n n=− + ⋅ [][][] 2d⎜⎟ 21 1TT−+ () − 1T⎝⎠ () d) x y For the analog system ∫ 2yyx+= ∫∫ -2 ∫ x[n] y[n] For the digital system from a) 2T− baae===−1; 1; 001 -2TD 2T− 1yne yn xn−−= [ ] [] e [] -1 x[n] (1+2T)y[n] For the digital system from b) T ; 1 2 ; 1bTa Ta==+ =− 00 1 12 1Tyn yn Txn+−−= [] [ ] [] D () 1 y[n-1] 9 For the digital system from c) 1 bb==1; 1; 01 (y[n]x[n] 21 T+ ) 21 ; 21aTaT=+ =− () () 01 21 21 1Tyn Tyn+ +− −= () [][] D D () 1xn xn=+ − [] [ ] -2(1-T) e) For ()sinxtt= and a sampling period of 0,1T= we have to find the amplitude of the output signal, in the analog and in the digital case sin ; 1 sin arg 1xt t yt H t H== + () () ()() {} () 11 HHs H== 1 ⇒=⇒ ω () ()() ω 22sjjj= ++ ω 11 HHj==1 ; arg 1 arg 2 arctg −+=− ()(){} {} 25 11 1⎛⎞ yt t⇒=sin arctg , amplitude − () ⎜⎟ 255⎝⎠ For the digital systems we have: 11nn⎛⎞⎧⎫ ⎛⎞ ⎛⎞ ⎛⎞ sin 0,1 sin sin argxn n yn H H== ⇒= + () [][] ⎬ ⎜⎟ ⎜⎟ ⎜⎟ ⎜⎟ dd⎨ 10 10 10 10⎝⎠ ⎝⎠ ⎝⎠ ⎩⎭ ⎝⎠ For the digital system from a) 111 () Hz=H⇒Ω= = () −− 21dT−−222djTTT−Ω − 1ez−⋅1cos sin1ejeez− Ω+ ⋅ Ω −⋅ 11 H⇒Ω= = () 242dTT−− 24212cosTTee−− +−Ω + Ω 1cos sinee−Ω () 11 1⎛⎞ H=≅ 2d⎜⎟ 0,2101− ⎝⎠ 0,4 0,2−− 1e− ee+− 12cos () 10 10 For the digital system from b) TTT Hz H=⇒ Ω= = ()() 1ddj−Ω − 12cos sin1212TjTzTe+Ω− Ω +⋅ +⋅ TT = = H⇒Ω () 222d Ω 14 4cos12cos 2sinTTTT++ +Ω+ Ω ()() 10,10,1⎛⎞ 0, 083H⇒≅ == d⎜⎟ 10 1, 21, 44⎝⎠ For the digital system from c) Ω − 1j− 11 11ze++ Hz H=⋅ ⇒ Ω=⋅ ()() − 1ddj−Ω 22111, 1 0, 9TTze+− − − () 22 1cos sin1122cos+Ω+ Ω () +Ω =⋅ =⋅ H⇒Ω () 22d 2 2 1, 21 0,81 1,98cos+ 1,1 0,9cos 0,9sin−Ω −Ω+ Ω ()() to 0,1 (homework) Set the value of Ω Problem 5 Consider a continuous-time system with the impulse response () and the aht transfer function () aHs The differential equation is: kkNN dy dx ab= kkkk∑∑ 00kkdt dt== We approximate this system using a discrete-time system using a method similar to the method of the differential equation approximation: 1xnxndx+− [][] ≅ tnTdt T= We define the difference of 0 ) ( order 0: [] xnxn∇= [] {} 11xnxn+− [][] ) ( order 1: [] xn∇= {} T 1kk− ()() 1 xnxn∇=∇∇ order k: [] [] {}{} {} s a) Let () Find the transfer function of the digital equivalent system Hs= 2a 1s+ () Hz ; () d b) What is the connection between ( )() aHs and Hz ? d 11 Solution 11xn xn+− [][] ) ( xn∇= [] {} T 11xn xn⎧⎫ [][] ) ( 2+− xn∇= ∇ [] {} ⎨⎬ T ⎩⎭ 111⎡⎤ )) (( 1xn xn=∇ +−∇ [][] {}{} ⎢⎥ T ⎣⎦ ⎤ 2111xnxnxnxn⎡ +− + +− [][][][] =− ⎢⎥ TT T ⎣⎦ 1 =+−++ 2 2 1xn xn xn⎡⎤ [][][] 2⎣⎦ T We notice that: k 1l kl +− 1xn Cxn l∇= () [][] {} kk∑ 0lT= a) We have to find the transfer function of the digital equivalent system ()Hz using a d method similar to the differential equation approximation bss 1 Hs== () 22a + + 21ss aasas++ 01 2 = = = = 0; 1; 1; 2; 1bbaaa⇒= 0101 2 The linear differential equation is: 2 dy t d y t dx t ()()() 2yt++ = () 2 dt dtdt The finite difference equation is: 12 1 ))) ((( 2ynynynxn+∇ +∇ =∇ [][][][] {}{}{} meaning that: +− 111ynyn xnxn+− [ ][][][] 2221yn yn yn yn⎡⎤ +++−++= [][][][] 2⎣⎦ TTT or: 2 21 221 1T n T n n n n n Tn Tnyyyy yyxx++−++−++=+− [][][][][][][][] () 12 Applying in both sides the Z transform, we obtain: 22 21 2TzTz zzzzz zTzzTzYYYYYXX+− + − += − () ( ) () () () () () () or: 22 22 21Yz T Tz T z z Xz Tz T⎡⎤+−+−+= − () ()( ) ⎣⎦ Tz T− Hz⇒= () 22d 21 21zTzTT+−+−+ () Hz ? b) What is the connection between ( ) ) aHs and ( d 1TzTz T− () − Hz== () 22d +− 11zT zT+− () () 1z− 11Tz− () T ⋅ T== 22 T ⎛⎞ ⎛⎞ 211zz−− 11T++ ⎜⎟ ⎜⎟ TT⎝⎠ ⎝⎠ 1z− T Hz= meaning that: () 2d 1z− ⎛⎞ 1+ ⎜⎟ T⎝⎠ We can see that: () = () 1zaHz Hs− ds= T 13